hep-ph/0303043 

CERN-TH/2003-051 
UMN-TH-2132/03 
TPI-MINN-03/07 



Supersymmetric Dark Matter in Light of WMAP 



John Ellis 1 , Keith A. Olive 2 , Yudi Santoso 2 and Vassilis C. Spanos 2 

1 TH Division, CERN, Geneva, Switzerland 
2 William I. Fine Theoretical Physics Institute, 
University of Minnesota, Minneapolis, MN 55455, USA 

Abstract 

We re-examine the parameter space of the constrained minimal supersymmetric exten- 
sion of the Standard Model (CMSSM), taking account of the restricted range of QcDMh 2 
consistent with the WMAP data. This provides a significantly reduced upper limit on the 
mass of the lightest supersymmetric particle LSP: m x & 500 GeV for tan (3 ^ 45 and \i > 0, 
or tan/? ^ 30 and fi < 0, thereby improving the prospects for measuring supersymmetry at 
the LHC, and increasing the likelihood that a 1-TeV linear e + e~ collider would be able to 
measure the properties of some supersymmetric particles. 
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1 Introduction 



The recent data from the WMAP satellite [1] confirm with greater accuracy the standard 
cosmological model, according to which the current energy density of the Universe is com- 
prised by about 73 % of dark energy and 27 % of matter, most of which is in the form 
of non-baryonic dark matter. The WMAP data further tell us that very little of this dark 
matter can be hot neutrino dark matter, and the reported re-ionization of the Universe when 
the redshift z ~ 20 is evidence against warm dark matter. WMAP quotes a total matter 
density Q m h 2 = 0.135j#oo9 and a baryon density Q h h 2 = 0.0224 ± 0.0009 [1], from which we 
infer the following 2-cr range for the density of cold dark matter: VLcdm^ 2 — 0-1126lo;oi8i- 
This range is consistent with that inferred from earlier observations [2, 3], but is significantly 
more precise. 

It has been appreciated for some time that the lightest supersymmetric particle (LSP) 
is a suitable candidate for this non-baryonic cold dark matter [4]. The LSP is stable in 
supersymmetric models where R parity is conserved, and its relic density falls naturally 
within the favoured range if it weighs less than ~ 1 TeV. This statement may be made more 
precise in the minimal supersymmetric extension of the Standard Model (MSSM), in which 
the LSP is expected to be the lightest neutralino x, particularly if the soft supersymmetry- 
breaking mass terms mi/ 2 ,rao are constrained to be universal at an input GUT scale: the 
constrained MSSM (CMSSM). Since Q x h 2 oc rn x n x , where n x is the relic LSP number density, 
and n x typically increases as the universal soft supersymmetry-breaking mass parameters 
are increased, one would expect the upper limit on m x to decrease when the upper limit 
on QcDAih 2 is decreased. Compared with taking this upper limit to be 0.3, as we and 
others [5, 6, 7] have done previously, taking Q C DMh 2 < 0.129 as suggested by the WMAP 
data - which is line with estimates from previous CMB determinations used in [8] - may 
therefore be expected to improve significantly the corresponding upper limit on m x . 

Such is indeed the case for tan/9 < 45 and fi > 0, or for tan/3 < 30 and /j, < 0, 
as we show below, where the largest values of m x are found in the x — t coannihilation 
region. We find for these cases that m 1( / 2 ^ 900 — 1200 GeV for QcDMh 2 < 0.129, whereas 
w-i/2 ~ 1400 — 1750 GeV would have been allowed for QcDAih 2 < 0.3. Correspondingly, the 
upper limit on the LSP mass becomes m x ^ 400 — 500 GeV, rather than m x ^ 600 — 700 GeV 
as found previously [9, 10]. This stronger upper limit improves the prospects for measuring 
supersymmetry at the LHC. Also, it would put sleptons within reach of a 1-TeV linear 
e + e~ linear collider, whereas previously a centre-of-mass energy above 1.2 TeV might have 
appeared necessary [11]. All the above remarks would also apply if other particles also 
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contribute to VLcdm- 

For any fixed value of tan/9 and sign of fx, only a narrow region of CMSSM parameter 
space would be allowed if 0.094 < Vt x h 2 < 0.129, as would be implied by WMAP if there 
are no other significant contributors to the density of cold dark matter. The previous 'bulk' 
regions of parameter space at small values of m Q and now become quite emaciated, and 
the previous coannihilation strips now become much narrower, as do the rapid-annihilation 
funnels that appear at larger tan/5 [12, 5]. However, unlike the coannihilation strips, the 
rapid-annihilation funnels still extend to very large values of m and mi/ 2 , the absolute 
upper limit on m x is much weaker for tan/? ^ 50 if fx > 0, or for tan j3 ^ 35 if fx < 1 . The 
narrowness of the preferred region implies that tan (3 could in principle be determined from 
measurements of mi/ 2 and m , as we discuss later. 

2 WMAP Constraint on the CMSSM Parameter Space 

Fig. 1 displays the allowed regions of the CMSSM parameter space for (a) tan/5 = 10, fx > 0, 
(b) tan/3 = 10, fx < 0, (c) tan/3 = 35,/i < 0, and (d) tan/3 = 50, ft > 0. We have taken 
m t — 175 GeV and A = in all of the results shown below. In each panel, we show the 
regions excluded by the LEP lower limits on rag, m x ± and rrih, as well as those ruled out by 
b — >• sj decay [14] as discussed in [15]. In panels (a) and (d) for fx > 0, we also display the 
regions favoured by the recent BNL measurement [16] of — 2 at the 2-cr level, relative to 
the calculation of the Standard Model based on e + e~ data at low energies [17]. 

Also shown in Fig. 1 are the 'old' regions where 0.1 < VL x h 2 < 0.3, and the 'new' 
regions where 0.094 < VL x h 2 < 0.129. We see immediately that (i) the cosmological regions 
are generally much narrower, and (ii) the 'bulk' regions at small m 1( / 2 and ra have almost 
disappeared, in particular when the laboratory constraints are imposed. Looking more closely 
at the coannihilation regions, we see that (iii) they are significantly truncated as well as 
becoming much narrower, since the reduced upper bound on Q x h 2 moves the tip where 
m x = nif to smaller m^. It is this effect that provides the reduced upper bound on m x 
advertized earlier. In panels (c) and (d), we see rapid-annihilation funnels that (iv) are also 
narrower and extend to lower rai/2 and ra than previously. They weaken significantly the 
upper bound on m x for tan (3 <^ 35 for fx < and tan (3 ^ 50 for fx > 0. 

We take this opportunity to comment on some calculational details concerning the rapid- 
annihilation funnels. Comparison with other studies of these regions [7, 18] has shown the 

1 Strictly speaking, there is also a filament of parameter space extending to large m x in the 'focus-point' 
region [13] at large m , to which we return later. 
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Figure 1: The (mi/2,mo) planes for (a) tan/3 = 10,fi > 0, (b) tan/3 = 10,/i < ; (c) 
tan/3 = 35,/x < 0, and (d) tan/3 = > 0. In each panel, the region allowed by the older 
cosmological constraint 0.1 < Q x h 2 < 0.3 has medium shading, and the region allowed by the 
newer cosmological constraint 0.094 < Q x h 2 < 0.129 has very dark shading. The disallowed 
region where m fl < m x has dark (red) shading. The regions excluded byb-^s^ have medium 
(green) shading, and those in panels (a,d) that are favoured by g^ — 2 at the 2-o level have 
medium (pink) shading. A dot-dashed line in panel (a) delineates the LEP constraint on the 
e mass and the contours m x ± = 104 GeV (rrih = 114 GeV) are shown as near-vertical black 
dashed (red dot-dashed) lines in panel (a) (each panel). 
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importance of treating correctly the running mass of the bottom quark, as we have done 
in previous works [5]. Also important is the correct treatment of annihilation rates across 
the convolution of two Boltzmann distributions when the cross section varies rapidly, as in 
the rapid-annihilation funnels, and we have taken the opportunity of this paper to improve 
our previous treatment. The results on the low-mi/ 2 sides of the the rapid- annihilation 
funnels are indistinguishable from those shown previously, apart from the narrowing due to 
the smaller allowed range of Vt x h 2 . However, there are more significant differences on the 
high-mx/2 sides of the the rapid- annihilation funnels, where our previous approximation was 
less adequate. This is most noticeable for the case tan/3 = 50, /x > shown in panel (d) of 
Fig. 1, where the two strips with relic density in the allowed range are all but merged. A 
second side of the rapid-annihilation funnel becomes distinctly visible when tan/3 > 51, but 
the gap between the two sides is much narrower than we found previously. This is also true 
for /i < 0, and is exemplified in Fig. lc for tan/3 = 35 where both sides of the funnel region 
are clearly distinct 2 . 

Before discussing further our results, we comment on the potential impact of a re- 
evaluation of m t . A recent re-analysis by the DO collaboration favours a central value some 
5 GeV higher than the central value m t = 175 GeV that we use [19]. If confirmed, this would 
shift the displayed contours of m h = 114 GeV to lower mi/ 2 , e.g., from ~ 300 GeV to 
m-i/2 — 235 GeV in the case tan/3 = 10, /x > 0. This would not affect the upper bound on 
m x that we quote later, but it would weaken the lower limit on m x in cases where this is 
provided by the LEP Higgs limit. 

The focus-point region would, however, be affected spectacularly by any such increase 
in m t , shifting to much larger mo- In the analysis of [13], for tan/3 = 10 and mi/ 2 = 300 
GeV, the focus-point is pushed up from m ~ 2200 GeV to m ~ 4200 GeV when m t is 
increased from 175 GeV to 180 GeV, and for tan/3 = 50, it is pushed up from m ^ 1800 
GeV to m ^ 3000 GeV. In our treatment of the CMSSM, we in fact find no focus-point 
region when m t = 180 GeV. In view of this instability in the fixed-point region, we do not 
include it in our subsequent analysis: our limits should be understood as not applying to 
this region, though we do note that it would also be further narrowed by the more restricted 
range of Vl x h 2 . 

We display in Fig. 2 the strips of the (mi/ 2 ,mo) plane allowed by the new cosmological 

constraint 0.094 < Vt x h 2 < 0.129 and the laboratory constraints listed above, for /x > and 

values of tan/3 from 5 to 55, in steps A(tan/3) = 5. We notice immediately that the strips 

2 Note that the irregularities seen in the cosmological regions are a result of the resolution used to produce 
the figures. 
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Figure 2: The strips display the regions of the (toi/2,TOo) plane that are compatible 
with 0.094 < Q x h 2 < 0.129 and the laboratory constraints for /i > and tan/3 = 
5,10,15,20,25,30,35,40,45,50,55. The parts of the strips compatible with — 2 at the 
2-a level have darker shading. 

are considerably narrower than the spacing between them, though any intermediate point in 

the (mi/2,mo) plane would be compatible with some intermediate value of tan/?. The right 

(left) ends of the strips correspond to the maximal (minimal) allowed values of m\/2 and 

hence m x 3 . The lower bounds on m 1( / 2 are due to the Higgs mass constraint for tan/3 < 23, 

but are determined by the b — > 57 constraint for higher values of tan (5. The upper bound 

on rrii/2 for tan/3 ^ 50 is clearly weaker, because of the rapid-annihilation regions. 

Also shown in Fig. 2 in darker shading are the restricted parts of the strips that are 

compatible with the BNL measurement of — 2 at the 2-a level, if low-energy e + e~ data 

are used to calculate the Standard Model contribution [17]. If this constraint is imposed, 

the range of mi/2 is much reduced for any fixed value of tan/3, and in particular the upper 

bound on ui\i% is significantly reduced, particularly for tan/3 50. However, there is in 

3 The droplets in the upper right of the figure are due to coannihilations when f is sitting on the Higgs 
pole. Here this occurs at tan/3 = 45. 
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general no change in the lower bound on mi/ 2 . 

3 Improved Upper Limit on the LSP Mass 

We now draw some conclusions from Fig. 2. Its implications for the allowed range of the 
LSP mass m x as a function of tan/3 are displayed in Fig. 3. As already mentioned, the 
upper limit is rather weak for tan/3 ^ 50 when /x > 0. However, for tan/3 < 40, we find the 
absolute upper bound 

m x < 500 GeV, (1) 

to be compared with the range up to ~ 650 GeV that we found with the old cosmological 
relic density constraint. Also shown in Figs. 2 and 3 is the strengthened upper bound on 
rai/2 and m x that would apply if one used the g^ — 2 constraint. We find 

m x < 370 GeV, (2) 

for all values of tan/3. Fig. 3(a) also shows the lower bound on m x as a function of tan/3, 
leading to 

m x > 108 GeV, (3) 

for all values of tan/3, with the minimum occurring around tan/3 = 23, when the b — > sj 
constraint begins to dominate over the Higgs mass constraint. As such, this lower limit 
depends on the calculation of rrih, for which we use the latest version of FeynHiggs [20]. 
This calculation has an estimated theoretical uncertainty ~ 2 GeV, and is very sensitive 
to m t . The lower bound (3) would become > 86 GeV (with the minimum occurring at 
tan/3 ~ 18) if we used m t ~ 180 GeV in FeynHiggs, or > 84 GeV (with the minimum 
occurring at tan/3 ~ 17) if we allowed for a 2 GeV reduction in the calculated value for the 
nominal value of m t . 

We do not show the plot corresponding to Fig. 2 for /i < 0, but we do show in Fig. 3(b) 
the corresponding lower and upper bounds on m x for tan/3 ^ 40. We note again that the 
upper bound would rise rapidly for larger tan /3, due to the appearance of a rapid-annihilation 
funnel analogous to that appearing for tan j3 ?Z 50 when fi > 0. In the fx < case, there is no 
possibility of compatibility with — 2 when the e + e~ data are used. We find the following 
range for tan/3 < 30: 

160 GeV <m x < 430 GeV (4) 
for ii < 0, with the lower bound being provided by b — >• 57 for tan j3 > 8. 
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Figure 3: The ranges of m x allowed by cosmology and other constraints, for (a) fx > and 
(b) fx < 0. Upper limits without (red solid line) and with (blue dashed line) the — 2 
constraint are shown for fx > 0: the lower limits are shown as black solid lines. Note the 
sharp increases in the upper limits for tan/3 ^ 50, /i > and tan/3 <i 35, fx < due to the 
rapid- annihilation funnels. Also shown as dotted lines are the e~L and x masses at the tips 
of the coannihilation tails. 
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We note that the upper and lower limits meet when tan (3 = 3.5(4.3) for fi > (<)0, 
implying that lower values of tan/? are not allowed within our analysis of the CMSSM. This 
lower bound on tan (3 is strengthened to 7 if the — 2 constraint is included. 

4 Implications for Supersymmetric Phenomenology 

The reduced upper limit on the sparticle mass scale improves the prospects for measuring 
supersymmetry at the LHC. Also, the fact that only a narrow strip in the (m^jWo) plane 
is allowed for each value of tan /? offers the possibility of determining tan (3 once mi/ 2 and m 
are known, e.g, from measurements at the LHC. We have discussed previously the sensitivity 
of Vt x h 2 to variations in the CMSSM parameters [21], and that analysis can be adapted to 
the present situation. Since the typical separation between strips with A (tan/?) = 5 is 
Am ~ 25 GeV and the width of a typical strip is Am ~ 5 GeV, it would in principle be 
possible to fix tan/3 with an accuracy A (tan/?) ~ 1 using measurements of m 1 / 2 and m 
alone for a fixed value of A (taken to be here). The required accuracy in mi/ 2 is not very 
demanding, since the strips are nearly horizontal, but mo would need to be determined with 
an accuracy Am ^ 5 GeV. It is interesting to compare with the accuracies in m , mi/2 and 
tan/? reported in [22] for the case m = 100 GeV, mi/ 2 = 300 GeV and tan/? = 2.1. This 
value of tan/? is not compatible with our analysis, but the expected accuracies A(mo) ~ 
10 GeV and A(m 1 / 2 ) = 20 GeV would already allow interesting crosschecks of the value of 
tan (3 extracted from a fit to the LHC data with that required by cosmology. 

The strengthened upper bound (1) on m x may also have important consequences for 
linear e + e~ collider physics [23]. At the tip of the coannihilation region, which corresponds 
to the upper bound in Fig. 3, we have m x = m fl , with the p,R and e# not much heavier. 
Therefore, a linear e + e~ collider with centre-of-mass energy 1 TeV would be able to produce 
these sleptons. This conclusion holds only if one restricts attention to the CMSSM, as 
studied here, ignores the focus-point region as also done here, and discards large values of 
tan /?. Moreover, we note that the left-handed sleptons are somewhat heavier at the tip of 
the CMSSM cosmological region: rrig ^ 700 GeV, as shown by the pale blue dotted lines 
in Fig. 3, and the lightest chargino has m x ± <; 800 GeV (green dotted lines). However, the 
strengthened upper limit on & x h 2 that has been provided by WMAP does strengthen the 
physics case for a TeV-scale linear e + e~ collider, compared with [11]. 
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5 Perspective 



In general, the WMAP constraint on Q x h 2 put supersymmetric phenomenology in a new 
perspective, essentially by reducing the dimensionality of the parameter space: one can now 
consider to be (almost) fixed in terms of the other parameters. The 'Snowmass lines' [24] 
now intersect the allowed cosmological region in just one (fuzzy) point each. The post-LEP 
benchmark points [25] have values of Q x h 2 that lie above the WMAP range. However, most 
of them can easily be adapted, in the 'bulk' and coannihilation regions simply by reducing 
itlq. An exception is benchmark point H, which was chosen at the tip of a coannihilation 
tail: WMAP would require this to be brought down to lower m^, which would make it 
easier to detect at the LHC or a future linear e + e~ linear collider. A more detailed update 
of the CMSSM benchmarks will be presented elsewhere. 
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